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1 Introduction 

> 

^ '. Recursion operators A± play an important role in the theory of the nonlinear evolution 

! equations (NLEE), integrable by the inverse scattering method. They have been con- 

structed and analyzed for a wide class of Lax operators L and appeared to generate not 
^ ' only the Lax representations, but also the hierarchy of NLEE's related to a given Lax 

operator L, their conservation laws and their hierarchy of Hamiltonian structures, see 
[31 EU [221 Is] and the numerous references therein. Such operators can be viewed also as 
the Lax L operator, taken in the adjoint representation of the relevant Lie algebra g. 
^ \ The construction of K± for Lax operators, whose explicit dependence on the spectral 

parameter A is comparatively simple (say, linear, or quadratic) has been done a long time 
ago [H [121 in Ej- Furthermore, the completeness property for the set of eigenfunctions of 
A± (the 'squared solutions' of L) is of paramount importance. The completeness of the 
'squared solutions' plays a fundamental role in proving that the ISM is, in fact, a nonlinear 
analogue of the Fourier transform, which allows one to linearize the NLEE. Using these 
relations one is able to derive all fundamental properties of the NLEE on a common basis. 

An important tool of extracting new NLEE from a known multi-component ones, is the 
the reduction group, introduced in [HI [151 [16]. It led to the discovery of the 2-dimensional 
Toda field theories [Ml [20] . Latest developments of the method led to the discovery of 
new automorphic Lie algebras and their classification [131 UHl US] • 

The problem of deriving recursion operators becomes more difficult when we impose 
additional reductions on L. If this additional reduction is compatible with L, being linear 
or quadratic in A, the construction of A is not a difficult task, see [H]- An alternative 
construction of A as formal recursion operator is given in [ini E] , see also the review papers 
[21 [21] . The effect of the Z,„-reduction is as follows: i) the relevant 'squared solutions' have 
analytic! ty properties in sectors of the complex A-plane closing angles tt/tt,; ii) the grading 
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of the Lie algebra q = ©^^qS'-'^^ is more involved and as a consequence the recursion 



operator is factorized into a product of n factors A = YYk=o ^k, and each of the factors 
Afc maps Afc: g^'^"^) ^ g(*^). 

The situation becomes more complicated if the additional reduction drastically changes 
the A-dependence of the Lax operator. Here we address this problem for one of the simplest 
nontrivial cases when the Lax operator L, due to an additional Z2-reduction changes its 
A-dependence from polynomial into rational one. The linear in A Lax operator 

LoV^ = z9^V^ + ALiV'(x,A) = 0, (LI) 

was shown to give rise [T0| |6] to the integrable system 

iut = {{1 - uu^)u^)^ , u^u = TL, (1.2) 

where u is (A^ — k) x k complex matrix and U is the unit matrix. The system (11. 2p is 
S{U{N — k) X U{k)) invariant, and, in this sense, is isotropic. In particular, if /c = 1, 
equation (II. 2p can be seen as a U{N — 1) invariant integrable system on CP^~^. The 
corresponding recursion operators, relevant for the fundamental properties of eq. (II. 2p . 
were derived in [7]. 

In [B] we derived new NLEE with additional Z2-reduction, which maps A — A^^. This 
reduction is among the simplest nontrivial ones, see [HI [T51 [ISl HSl [IS]- Due to it, the 
relevant Lax operator acquires a rational dependence of A: 

Lij = id^ij+(xLi + \L_i]4j = 0, (1.3) 



;i-4) 



A 

and the simplest NLEE take the form 

iut = - {u{u*u,j, + v*v^))^ + %vv*u 
ivt = Vxx — {y{u*Ux + v*Vx))x — 8uu*v 

where u and v are functions of x and t subject to the condition: 

\u'^\ + \v^\ = l (1.5) 

i.e. the vector with components u and v sweeps a 3-dimensional sphere in M^. System 
(11.41) can also be seen as an anisotropic deformation of (11.21) with k = 1, N = 3. 

Our main aim in the present paper is to present two ways of deriving the recursion 
operators related to the Lax operators with rational dependence on A. It is a natural 
extensions of our previous results [6l [7] . In the next Section 2 we give preliminaries con- 
cerning the spectral properties of L i \1.3\\ . In Section 3 we derive the recursion operator 
using the Giirses-Karasu-Sokolov (GKS) method. In Section 4 we use the Wronskian re- 
lations to determine the 'squared solutions of L. Next using the gauge covariant approach 
[H [8] we introduce the recursion operator as the one for which the 'squared solutions' 
are eigenfunctions. In the last Section 5 we give some conclusions, while the Appendix 
contains some details from our calculations. 



2 Preliminaries 

In this section we first formulate the basic results from our previous paper |6]. There 
we have shown that eq. (11.41) . which is naturally related to the symmetric space [11] of 
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A. Ill-type SU{3)/ S{U{1) x U{2)), allows a Lax representation and can be solved by the 
inverse scattering method. The Lax operator L is given by eq. (11. 3p and the equation 
itself is the compatibility condition between L and the linear operator 

Ai; := tdti; + + AAi + + X^A^ + ^A_2^ ^ = ^f{X), (2.1) 

It is well know the above mentioned symmetric space is constructed by using Cartan's 
involutive automorphism [H], which induces a Z2 grading in the underlying Lie algebra 

g = g(0)©g(l), Ji = diag (1, -1, -1) 

= {F e 0: J^YJ^ = Y}, g« = {X e g: -hX-h = -X}. 

Our Lax representation is such that L±i, A±i G and Aq, A±2 £ 0*^, for more details see 
the Appendix. In particular, for Li and Ai we write 



Li = n* , Ai = a* I . (2.3) 
\v* Qq) \ 6* 

which explicitly involves the first two reduction conditions: 

4i = L±i, 4i = ^±i, 42 = ^±2. (2.4) 

We impose also a third reduction of the form 

J2L1J2 = L_i, J2A1J2 = A_i, J2A2J2 = A_2 (2-5) 

where J2 = diag(l, —1, 1). At this point, we impose one additional requirement 

that is the vector {u^v) lives in a 3-dimensional sphere in R^. 

The spectral theory of the Lax operator L (II. 3p substantially depends on the boundary 
conditions. We have two natural ways to fix the boundary configurations, namely: 

1. \mi^-,±oou{x,t) = 1, lim^_^±^v{x,t) = 

2. lim^^±oou{x,t) = 0, lim^^±oov{x,t) = 1 

The Jost solutions of L are determined by using the asymptotic potentials U±^as{X) = 
lima;_j.±oo(AI/i+A~^L_i) which can be diagonalised for both choices of boundary conditions 

1 1 

1) 9o'U±,U^)9o = Ja(A), g, = ^\ 1 -1 




-V2 
J,(A) = (A -A-i)diag (1,0,-1) 

^ /I -1 

2) go'U±,M)go = MX), ^0 = ^ 

\1 1 

Jb(A) = (A + A-^)diag (1,0,-1). 



(2.6) 




The continuous part of the spectrum of L for the case a) fills in the real axis M while for 
the case b) it is the union of the real axis and the unit circle, see FiglH 

In [6] we have constructed the fundamental analytic solutions (FAS) of L for both 
cases. Skipping the details, we note only that these FAS can be viewed as solutions of 
Riemann-Hilbert problems (RHP): 

The Jost solutions of L and the related scattering matrix T(A) are defined as follows: 

hm ^±(a:,A)e-^^W^^7o'i = Il, 

r(A) = ^;i(x,A)^_(x,A) 

For the case a) the fundamental analytic solutions ^ind x~ ^'^^ solutions of the following 
RHP 

X+(x,A) = x-(x,A)G(A), G{X) = [S-{\)r^S+{\), A G R. (2.8) 
They obey the symmetry relations imposed by the reductions 

(x+)^(x,A*) = x-(x,A) 
JiX+(x,-A)Ji = x-(x,A) (2.9) 

In the case b) we have 4 fundamental analytic solutions x^^^ with analyticity regions 
which satisfy a RHP on the contour shown on Figure 1: 

X^'\xA)=X^^\xA)G{\), AG [-1,1] 

X^^{x, A) = x^^\x, A)G(A), A G (-oo, -1] U [1, oo) 

X^^\xA)=X^'\xA)G{\), A = e^^ ^e(0,7r) 

X^'\xA)=X^^\xA)G{\), \ = e'\ (^G(7r,27r) 

where the sewing function is given by 

G(A) = {S~{\)r'S\\). 

The FAS obey the following symmetry relations: 

[x(i)(x,A*)]^ = \^^\x,X)Y\ [x^^)(^,A*)]^ = [x(=')(x,A)]-^ 

Ja^^)(x,-A)Ji =x('n^,A), hx^'^\x-X)h = ^\x,X) (2.11) 

J2X^i)(x,l/A)J2 = X^'na^,A), J2x(2)(x,l/A)J2 = X^'H^,A). 



(2.10) 
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3 Recursion operators 



This section is dedicated to the construction of recursion operator for the NLEE f ll.4p . 
For this to be done we are applying the method proposed by Giirses, Karasu and Sokolov 
(GKS) [lOj. The recursion operator is viewed here a pseudo-differential operator to map 
a symmetry of a NLEE into another symmetry. The underlying idea of the GKS method 
is to interrelate two adjacent flows A and A with evolution variables r and t respectively 
through a equality in the form [101 [23] 

V = {X"^ + X-^)V + B. (3.1) 

The matrix-valued functions V and V are involved in the Lax operators A and A as 
follows 

A = idr + V{x,t,X), (3.2) 
A = idt + V{x,t,X) (3.3) 

All quantities above must be invariant under the action of the Z2 reductions (12. 4p and 
(12.51) . This explains the choice of the coefficient before V in (13. ip which is invariant under 
A — !■ —A and A — A~^. The operator B is chosen to be a rational function of A 

B = Bo + XBi + jB_, + X^B2 + ^B_2, B_k = J2BkJ2, k = 1,2 (3.4) 
The matrices involved in B are hermitian ones and have a block structure given by 

--(c-oi' --(o^)' -K;-)' 

in a complete analogy with the second Lax operator (12. ip . 

At the begining of our consideration, let us remark that the zero curvature condition 
of the Lax operators L and A could be written in the following manner: 

tLt = mV], (3.6) 

where Lt = XLi t + X^^L^i t- After substituting (13. ip in the analog of (13. 6p . when the 
evolution paramenter t is replaced by r, we obtain the following basic equation: 

tLr = i{X^ + X-^)Lt + [L,B]. (3.7) 

Then the recursion operator TZ acts as follows: 

Since all matrices are traceless, we have 

d= -tiD = -{a + S). 
By comparing coefficients before different powers of A, one splits equation (13. 7p into 

iLi^t + [Li, B2] = lut + (D* - d)u = 0, (3.9) 
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iB2,. + [Li,Bi] = 



iD^ + u*c^ - c*u^ = 





id.j. + c* — c'^u* 



iLi,, = iL^i^t + iBi,,, + [Li, Bo] + B^] 
iBo,.+ + =0. 



(3.10) 

(3.11) 
(3.12) 



The rest of relations could be obtained from those above, by multiplying each quantity 
by J2 on its left and on its right hand side. 

The system (13. 9p is linear for the matrix elements of D. A solution to (13. 9p is given 

by 



a = i{uu* + v*Vt) — d{2\v\^ — |mP) (3 = i{vul — u*Vt) + 3du*v 
S = -i{uul + v*vt) - rf(2|up - \v\^). 

The solvability condition for the system of linear equations (I3.13P is: 



ttxl'yp + ^zkl^ — (3xUV* — (3*u*v = 



One can check that 



-v{5x - a^) - 2u(3x 



K 



(3.13) 



(3.14) 



(3.15) 



where k is real, but otherwise arbitrary, is the general solution of fl3.13p . In what follows 
we shall set k = without any loss of generality. It could be written in the following 
matrix form 



s 
c* 



A 



/ /3 \ 

a— 5 

2 ^ 
, a+o 
\ 2 / 



A := i — 

dx 



/ -V -u 0\ 

-u V 
V* u* 
\ u* -V* J 



dx 



(3.16) 



Next step is to find the function d, involved in the expressions for a, P etc. For this to 
be done, we make use of the condition (I3.14p which leads to a linear differential equation 
for d: 



- 24. + i{u*ut + v*vt)x - i [{uut + v*vt) 



\v\ — \u\ 



- {uvl - v*ut){u*v)x + {u*vt - v*ut){uv*)x] 
After some simple manipulations one obtains 



0. 



-1 ._ 



d = --{uu* + vv^) - -d^^ {utul + vtvl - u*Ux - v*Vx) , 



(3.17) 



(3.18) 



where d'^ 

To calculate TZ, one needs the diagonal matrix Bq. It directly follows from fl3.12p that 
the following equation holds true 



uc* — u*c + v*s — vs* 
iBo^x = 2 I u*c-uc* 

vs* - v*s 



(3.19) 
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Taking into account formulas (13.101) . one deduces that 



Bn 



'a- 6 0' 
-a 
6 



Furthermore, taking into account the structure of (13. lip , one can easily see that TZ can 
be split in the following form: 

n = C + no + V. (3.20) 

The first term C originates from the first term in (13.111) . and it is a constant matrix 
diag (—1, 1, —1,1). The third operator V is obtained from the two commutators in (13. lip . 
It turns out that it splits into a local term 



loc 



/ 8|m|2-1 


-4(M*)2(3|t;|2-l) 
\ -12u*v*\v\'^ 



-4^2(3 |y I 



-12uv\v\'^ 

12 1 



12nt;|n|2 \ 



\u\ 







4{v*f{3\u\ 



and a nonlocal one 



noni 



f u{l-3\v\^) \ 
v{3\u\^-l) 
u* (3|tf -1) 

\v* {1 - 3\u\'^) J 



9x [«,<,-Mx,-^'x)- 



Finally the operator TZo = ^(-Bioc + -Bnoni) sphts into a local and a nonlocal parts as 
follows: 



/ 



loc 



4 







-Au* 2v{3\v\^-l) -6u*v^ \ 



2v*(l-3\vn 6u(v 



*\2 







4u 



-4u* 
-u* 



3n(|tf - Inp) -t;(l + 6|n|2) 



nonl 



3t 



/ 2u*v \ 
2uv* 

V -1/3 / 











(3.21) 



(3.22) 



In fact, TZq represents the recursion operator in the case of polynomial bundle Lax pair 

L := idx + XLi, 
A := tdt + AAi + 



4 Recursion operators and 'squared solutions' 

This approach is based on the Wronskian relations, mapping the potential of the Lax 
operator L onto the scattering data, which allow one to introduce the 'squared solutions' 
for L (see, e.g. [S] and the references therein). 
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The Wronskian relations, derived in [6], take the form: 

POO 

{iX~^Jox{x,\)-iJo,Eo,)\'^^_^= / ([Li, Jo],$i(x, A)) 



x=~oo 
oo 



i {x ^Sx{x, A), |^_^ = - / dx {6Li, <l>i(x, A)) , 



(4.1) 



where we have introduced the squared solutions: 

<^i{x, A) = Ae«(x, A) + A-Vo(ea)(a;, A), (4.2) 

and 

e„(x, A) = xEaX'^ix, A). (4.3) 

where is one of the Cartan-Weyl generators. 

In this approach, see [HIH], one picks up a certain 'projections' of the 'squared solutions' 
(see K^'^{x, A) in eq. ( 14.12^ below) which contribute to the right hand sides of eq. ( 14. ip . 
Then the recursion operators can be introduced as the ones, that have K^'^{x,X) as 
eigenf unctions. 

Let us introduce 

Ji = diag (1, -1, -1), K, = diag (1, 0, -1), Ko = diag (1,-1, 1) (4.4) 

and 

^) = ^) + ^"Vo(e^(x, A)), (4.5) 

which satisfy the equations: 

+ [^1' ^)] + [^"1' J^' ^)] = 

MK,ai^^ A)) = ^t,Jx, A), M^o-J^, A)) = $o1«(^, A) (4.6) 

In addition, we have: 

(A + A-^)<f A) = $^+1^ Jx, A) + ^tij^. A) (4.7) 
In what follows we will use only the equations for ^^.^{x, A) and ^f.^{x, A): 



. , :il + v.o)[i^i,<^'t(^'A)] = 

- [^1 - ^-1, <J>t(^' A)] = -(A + X-')[L„ A)] 
Next, we insert the splitting: 

$± Jx, A) = H^^^ix, A) + K^Jx, A), A) = H^Jx, A) + A) 

H^Jx, A) = ^o(i^o1a(^, A)), K^Jx, A) = v:^o(i^o1a(^, A)) 



(4. 



(4.9) 
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and obtain 



z-^ + (U + vo)[Li, H^Jx, A)] = 



dx 

[L^ - K^Jx, A)] = -(A + \-')[Lu K^Jx, A)] 



dx 

■• -[Li-L„i,/7o%(x,A)] = -(A + A-^)[Li,i/± (x,A)] 



(4.10) 



(4.11) 



dx 

The integration the first two of the above equations gives: 

HU^, A) = /io±o;« + ^(1 + ^^0)5,7^^1, i^it^(x, A)] 
Kti^, A) = A:o±o^„ + ^(11 + ^o)d;'[L„ Hf^{x, A)] 

where h^^.^^ and /c^.„ are matrix-valued constants. Thus we expressed K^.^ and H^^^ in 
terms of K^.'^ and -fff':^"'" where 

H^.:^^{x, A) = iJi%(x, A) - ^2) 

1 (4.12) 
Kt^{xA) = K^Jx,X) - -Li(i^J„,Li) 

where L2 = Li — 2/311. Let us now derive the equations for K'I:'J~ and H^. '^. To do this, 
we insert eq. (14.121 and insert it into the last two equations of (I4.10p . Thus we get: 

^ y^2,x(i^i;„, L2) + —L2 ^ [Li - Kl,^[x, A)] 

= -(A + A-i)[Li,AtJx,A)] ^^^^^ 

2 2 ' f ' " ^-i.^o%(^, A)] 

= -(A + A-^)[Li,i7±Jx,A)]. 

Next, we separate the left hand sides of eq. fl4.13p into terms commuting with Li and 
L2 and terms 'orthogonal' to them. In order to calculate the coefficients {H^^, L2) and 
{Kf:^,Li), we multiply both sides of the first (resp. the second) of the eqs. (I4.13P by L2 
(resp. by Li) and take the trace. The result is: 

■ / ^H^'a J \ , d{Hf L2) T^± / \\} r \ n 

" 'd^' 7 ' — ' ~ " 

Integrating fl4.14p we get: 

{H^^^, L2) = h^,^^ - d-' / L2)- id-\[L, - L_i, i^o^Jx, A)], L2) 



(4.14) 



(4.15) 



where k^^ and h^^ are matrix-valued integration constants. 
Skipping the details we obtain: 



-ad^j^([Li -L„i,A;oo]) 

where /loo and /cqo are arbitrary constants which we shall set to be equal to zero. The 
operators Aj, j = 1, 2 are defined as follows: 



-1 ^Kt, 



dx 



-zad 



-1 9H^. 



(9x 

-1 f/9Ht^ 



;i-^o) 



5x 



Thus, we obtain: 



A2Aii^&^ = (A + X-^Kt"- 
KrK2Ht'^ = (A + X-^H^.^^ 



(4.17) 
(4.18) 



5 Conclusions 

We derived the recursion operators A for the system (11.41) related to the symmetric space 
5'f/(3)/5'(f/(l) X U{2)) with an additional Z2-reduction. Our first derivation of A is based 
on the GKS method [10]. 

Another way to construct A consists in using the Wronskian relations. They allow us 
to determine the 'squared solutions' of L which are treated as eigenfunctions of A. 

Our results could be extended for operators L related to a generic symmetric space 
of the type A. Ill = SU{n + k)/S{U{n) x U{k)) as well as to other types of symmetric 
spaces. 
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A Some intermediate results 



In this appendix we will present some intermediate results, used in section 4. The invo- 
lutions J induce a Z2 grading in the Lie algebra q as follows: 



= 5^°^ 

^(0) ^^Y:Y = JYJ~^}, 0(1) = {X:X = -JXJ-^} 




ab 

Y = \ A;m|, X=|a*00 

6* 



(A.l) 



The inner automorphism ifo staying in the formula for the "extended squared solutions" 
(I4.9p can be represented in the form: 





Next we take into account that ipo{Ho) = Hq and (Po{Kq) = Kq and therefore: 

{[L^ - L^,,K^Jx,X)],L2) = 2\u\\4,v* - /^J^O 
([Li-L„i,i7ot(x,A)],Li)=0 



where 



Then 



/ dH^'^ \ 



00,a) 



(A.2) 



(A.3) 



i^o":. = I , Ht=\ fco±„ I , (A.4) 
JaO / V < 

Introduce the notations 



f-kta-nt^ 

Kt = I ^£ o h =\ I (A.5) 
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